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Механ!зм прискорення обчислень в метод! Литтла 
для розв’язку задач класу ком1вояжера 


Поиск решений задач класса коммивояжера в бинарной схеме ветвлений метода ветвей и границ можно 
значительно ускорить за счет обращении к быстрому алгоритму решения одного из вариантов задачи 
о назначениях (ЗН), применяемой для вычисления нижних оценок стоимости гамильтоновых маршрутов. 


Оптимальное решение ЗН для полученной матрицы можно найти за время О [8 . Предложен алгоритм по 


схеме ветвей и границ, использующий быстрый алгоритм решения ЗН в качестве нижней оценки стоимости 
решения. 
Ключевые слова: задача коммивояжера, точный метод, модификация метода Литтла. 
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Пошук розв’язюв задач класу компвояжера в бнарнйй схем! розгалуження метода плок та меж можна 
суттево прискорити за рахунок звернення до швидкого алгоритму розв’язку одного з варантйв задач! про 
призначення, яка використовуеться для обчислення нижн!х ошнок вартост! гампльтонових маршрут. 


Оптимальний розв’язок задач1 про призначення для отримано! матриц! можна знайти за час 0(»?). 


Запропоновано алгоритм за схемою плок та меж, який використовуе швидкий алгоритм розв’язку задач! 
про призначення в якост! нижньо! ощнки вартост! розв’язку. 
Ключов! слова: задача комвояжера, точний метод, модифлкащя метода Литтла. 


Введение 


Рассмотрим следующую формулировку задачи о назначениях (ЗН). Задана матрица 


хх р . + . 5 
стоимостей Ве: ‚ где с; = при {= Л ис; Е Ху или с, =® при 1], 5] =Ьп, 
п 
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Го — множество неотрицательных целых чисел. Требуется найти оптимальную пере- 
п 
становку с = (с [],< [2],...,6 ["]) стоимостью <) =тт Ус. где п=(п[],л[2],... п[]) = 
и 
перестановка и столбцов матрицы, с. и, #=Ьл. ЗН может не иметь решения если С 
содержит элементы с, =, ГЕ 
В статье [1] предложена процедура 5/1’, позволяющая по перестановке © для 


матрицы С, найти за время, пропорциональное о (”°) оптимальное решение о,’ 


для матрицы С 


и полученной из С, заменой значения с, на бесконечность. 


Статья посвящена модификации метода Литтла, использующей ускоренную про- 
цедуру вычисления нижней границы стоимости оптимальных решений для ключе- 
вых задач класса коммивояжера. 


Ускоренная процедура вычисления 
нижней границы стоимости 


Рассмотрим следующий вариант ЗН: известно решение о = (с [1], [2] в №] для 


матрицы стоимостей С =[с, | ‚ по заданной перестановке с надо найти решение ЗН 
п 


сь(х, у], содержащее элемент с,,, и не содержащее элемент сур, Су, С (с бо] 


51. 


Двудольный граф (Х,У,0) в котором элемент с, #® матрицы С является 


ху ? 


ребром [Ь Л, Г= ], содержит решение ЗН о для С ‚, как совершенное паросочетание с 


минимальным суммарным весом ребер. Решение о, [о у) для матрицы Со, [1% у), 
если оно существует, также представлено, как совершенное паросочетание в дву- 
дольном графе (Хх ВИ '), ПИ {0 ыы) [а, ь]} | Йа — множество всех ребер, инци- 
дентных вершинам х и у. Задача нахождения перестановки б/у (х, у) решается за время 


0 ("*) методом поиска кратчайшего увеличиваю щего пути из вершины а в вершину р 


в двудольном графе (0) Способ построения перестановки б,‚,, не вклю- 


чающей элемент (х, у) о и перестановки б | у] содержащей (х, у) О ДЛЯ ВЫЧИС- 


ления нижних границ в схеме ветвления алгоритма Литтла является базовой компо- 
нентой метода, обеспечиваю щего поиск точных решений задач класса коммивояжера. 
Этот класс образуют задача коммивояжера (ЗК), общая ЗК (ОЗК) и гамильтонова ЗК 
(ГЗК) [2]. 


Пусть С = [8 — матрица стоимостей, где с; = О (о в) в ОЗК, с. =4; вЗКи 
П 
ГЗК, 4 = (а[1],4[2],...9["]), © (&) - соответственно обход и его стоимость в любой 


из этих задач. 
Нижняя граница ф, в корне дерева перебора равна минимуму С (с ) целевой 


функции ЗН для некоторой матрицы С. Ее можно приблизить к стоимости оф (8°) 
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* 
оптимального обхода е следующим образом. Вначале матрица С приводится по 


строкам, а полученный результат по столбцам. Затем матрица С приводится по 
столбцам, а полученный результат — по строкам. Из них выбирается наименьшая 


сумма пу = тт{ Ну, Но} (в общем случае и! * и»), и соответствующая ей матрица Су, 
по которой находится решение ЗН о. За нижнюю границу стоимости любого обхо- 
да, порождаемого матрицей С, принимается величина 

Фо = Но + Со (во), но = т {ну, №} (р 


Временные затраты на вычисление фу оцениваются величиной 0("). Для про- 


стоты положим в корне дерева перебора С = Су, в =су, и = Но, Фи =С (с ) р 


Наибольшей неопределенностью, замедляющей процесс вычислений в методах 
типа ветвей и границ, характеризуется этап выбора наиболее подходящего для включения 
в оптимальное решение элемента из исходных данных, причем число кандидатов не 
меньше порядка п матрицы С. Покажем, ЧТО В рассматриваемой модификации ме- 
тода Литтла предпочтительный элемент, выбирается только из и элементов решения 
ЗН. 


р. 


терминах задач класса коммивояжера обходом. Чем меньше К, тем лучше структур- 
ные свойства ЗН, применяемой в качестве релаксации для построения замкнутых 
маршрутов в графовых моделях. 


Рассмотрим взвешенный граф (у, [#1 ) с множеством вершин У и множеством 


Перестановку с в корне дерева перебора представим цикловым разложением 


у=1 К } ,К< | и/ 2 | . При К =1 о -— циклическая перестановка, называемая в 


> + 
дуг И ‚ в котором пара вершин (, у у) является дугой весом с,„ © К , если в мат- 
рице С с,=®. Цикловое разложение перестановки с разбивает множество И на 
подмножества И, содержащее все дуги циклового разложения, и И,, к которому 
отнесем все остальные дуги. Если удаляется дуга (х, У) #0, то, очевидно, о, =, 
С (©) © (с я Добавление этой дуги к И; либо приводит к образованию переста- 
нОвки © (х, У), С (с (5 У)) >С (с ), либо вообще исключает ее построение (рис. 1, в). 
Поскольку С (.„)< С (с (с, у) ДЛЯ [% у) И; , то из двух активных вершин (<, у)°) 


И ((. у}, порожденных корнем ©, вершиной ветвления должна быть (+, >}. 


Таким образом, корню © и вершине ((.. ›)°) соответствует одна перестановка о, 


после построения которой была выбрана дуга [% у) О ,. Поэтому дугу (а. у), ИНИ - 


циирую щую ветвление, следует выбирать из подмножества (|. 


Рассмотрим, как в предлагаемой модификации метода Литтла строится дерево пере- 
бора И ВЫЧИСЛЯЮ ТСЯ НИЖНИЕ Границы стоимости решений 8 в его вершинах. Если ДЛЯ 


матрицы стоимостей ГЗК ЗН не имеет решения, то ГЗК также не разрешима. В общем 
случае процесс ветвления начинается с построения перестановки с, цикловое разложение 


которой образуют к (с ) контуров, К (с )> 2, вычисления С (с ) и нижней границы 


5 (+*)- С (с )+ и в корне дерева перебора. 
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= 


6) в) 
Рисунок 1 —а), 6) множество циклов, построенное путем добавления к цикловому 
разложению перестановки ©’ дуги (х, у) е Ио ; в) добавление дуги (х, у), 
не содержащейся в цикловом разложении перестановки с’, исключает построение 
вообще какой-либо перестановки. Добавленная дуга (х, у) изображена утолщенной 
линией, тонкими линиями представлены дуги, порождающие множество циклов из дуги 
(х, у) и циклового разложения перестановки с’; пунктирными линиями обозначены 


дуги, исключенные из циклового разложения перестановки о 


При исключении из решения ЗН о элемента (х, У), инициализирующего ветв- 
ление, степень его влияния на стоимость оптимального обхода определяется оценкой 
Г(х, у). Чтобы определить (х, у), каждую дугу (К, /), удаленную из циклового раз- 
ложения перестановки с оценим величиной 

у (К, 1) = тт су +тпт с; —сц (2) 
Е ЕК 
Решение ГЗК о ‚ матрица стоимостей которой содержит с; =®, 1* ], может 


включать элементы (К, 1 ) с оценками у (К, | ) = о. Такие же элементы могут быть в 
решениях ЗН, полученных при вершинах ветвления в процессе построения обхода 
ты ЗК или оптимального маршрута ОЗК. Очевидно, если у (6; 1) = со ‚ то обход, соответ- 
ствующий решению ЗН, включает (К, 1 ) ‚а Ру У) == оо. Нетрудно заметить, что не 


существует обхода, для которого построено решение ЗН, если в цикловом разложе- 
нии этого решения все дуги какого-либо из контуров имеют оценку с или число г 


дуг (,1) с оценкой со больше п-К. 
Список элементов с оценкой © в решении ЗН и элемент (х, У) с оценкой 
Г(х, У), если он существует, находится процедурой ГАМС . 


Вход процедуры представлен матрицей стоимостей С = [5 | с элементами 
п 


ООО Е с. =Ю, или с, =%, #=], решением ЗН с = ((к,1) |1 <", 
[= 11, 2,.... п ), си * 90 для матрицы С ‚ числом контуров К (с) в цикловом разложе- 
нии перестановки си пустым списком К (с )= ((, Г) | у (К, Г) = о), РЕ [Е (с ). 

В процессе поиска элемента (х, у), инициирующего ветвление, процедура 


возвращает список в (с) элементов, не содержащих [8 у), формирует цикловое 
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разложение перестановки о на контуры уз \ = 1, К ‚ определяет их мощности и 
выбирает наименьшую среди них Кит . Искомый элемент [8 У) существует и при- 


нимает максимальное действительное значение оценки, определяемое из (2), если 
г<К‚ а или г<п- К. Процедура ГАМС включает такие действия 


Берт 
К(с)=0; г=0; М =0; 
К:=1 
уе д < п 4о 
Бет 
с(в):= ши (с, | 1=1}; 
с (Г) := ши (си |1 К); 
у (&,Г):=с(#)+с(1)-сы; 
Из (К,1)=0 Ше 
Берт 
В (с ):= К (в) (Е, Г); г :=г+1; 
Ис(К)=х Шеп 
юга ]е {1,2,....1-1,1+1,... п} @0 
сн :=®; 
Ис(Г)=х Шев 
Гога Ге {1,2,...К-ГК+1,.... п} 60 
бр 
еп д 
Шу (,Г)>М Шеп М :=7 (&,1); 
К:= +1; 
еп д 


Иг>п-К Фей М =; 


преобразовывать перестановку о в цикловое разложение {2 и = 1 к} 


и найти Ки = ШШ {2 м - К}; 
И Кии <Г Фе 
регт 
и :=1; 
УПВНе и < К 40 
Берт 
2, =7,-В (6); 
НИ 2, =© Ше М :=©; ее и:=и+[ 
ета 
ета 
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Ё М =о Фев 
Безт 


(х, у) := (1); 


сохранить все элементы и их оценки в контуре й. а? содержащем 


(% У) 


епа 
епа 


2 
Процедура ГАМС. за время О(п преобразует входную матрицу С в матрицу 
С’ того же порядка, что и С, и в зависимости от мощности списка К(с) устанав- 
ливает, является ли элемент (К, [ ) с максимальным действительным значением М, 


определяемым из (2), искомым элементом ( р У) с оценкой 


Г (+, у) = тах (= (Е еб -К (о )}. 


Если входная матрица содержит элемент [% у), инициирующий ВЕТВЛЕНИе, ТО 
вместе с матрицей С’ и списком в (с) процедура запоминает все элементы подмно- 


жества 7,, которому принадлежит (х, у). Вершина ветвления, представленная корнем 
дерева перебора, порождает вершины (<. у}, ((х. у)°) и список 0 (с } элементов мат- 


рицы С’, изменяю щих свои действительные значения на ®. В корне дерева 0 (о )= О. 


Найдем в вершине ((. у) нижнюю границу о (5. у) стоимости любого 
обхода, не содержащего элемент (х, у)ес -А (с). Очевидно, искомая граница оп- 
ределяется в результате преобразования решения ЗН с в решение ЗН ©, для ма- 

| У ! у + 
трицы С’, полученной из С’ заменой с, Е № на ®. Выполним процедуру нахождения 


о. 5/1 . Тогда 


ху 


, (>) С' (6 ей ру если построена перестановка бу, 
о |= ” | 


и Иначе, 


Ясно, ЧТо вершина ((.. ›}. | НЕ ПОДЛЕЖИТ ВЕТВЛЕНИЮ , ЕСЛИ ° (2, у) |= ® ИЛИ 


Построенная перестановка б , циклическая. Если перестановка 0, не циклическая, 


то верш ине (с. у} | поставим в соответствие с, ф ((-. у} ( (о р Число кон - 
туров К (с и в цикловом разложении 6, и списки в(с„)=К(о), 0(с„)= 
= О(с ) №) ( ра У) ‚ восстанавливающие по исходной матрице матрицу С Г 

Определим в вершине ((х, У) с) оценку снизу Ф [( У) с) стоимости любого об- 


хода, содержащего элемент (1% У) и все элементы списка К(с). Такой оценкой 
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является стоимость перестановки с [% У) ‚ включающей (х, У) и К(с ). Перестановка 
[е [% У) может совпадать или не совпадать с о. Для того чтобы выяснить, когда 


б [8 У) =о и как найти с [8 У) =о , обратимся к контуру 7.,, которому принад- 


а 2 
лежит дуга (и, У) ‚ соответствующая элементу, инициирующему ветвление. [Е й =Па. 
Обозначим и, число дуг контура, получивших при поиске дуги (х, У) оценку оо. 


Каждой такой дуге соответствует элемент из списка К(с ) . Очевидно, при п, = 
= [Е 4|-1 не существует обхода, который содержит [5 У), т.е. ф ((х, У) с) (рис. 2, а). 

Пусть и, < [2.|--2. Тогда о (х,у)=о, Ф((х, У) °) =С'(с<)+ы, С’(с)=С(с), 
К(с [1% У)) = К (с у (<>, У). Список О(с (% У)) формируется за время (п, добав- 
лением к 0(с) элементов, которые при изменении в матрице С’ своих действитель- 
ных значений на с устраняют контур с дугой [0% У) и максимально возможным числом 
дуг из А (с ) . Рассмотрим все случаи исключения подконтуров. 

| О(с [7% У)) = О(с ) (у, х) ‚ если дуга (х, У) не инцидентна в контуре 7, дугам 
с оценками, равными ©, и с, #0 в матрице С 1 (рис. 2, 6). Положивв С’с Суд =, 


2 > 
получим матрицу С для поиска в неи очередного элемента, инициирующего ветв- 
ление. 


2. Если контур 24 содержит цепочку дуг Е ХУ, У=У,...У>5У,, У) с 


оценками '/ нь х]* со, У (У Уьн) 700, У(у-1ьу,) =00, г=2,р,и сх = со (рис.2, в), 


то в матрице С' положим с, - =50. В результате получим матрицу С *и список 
р 


О (с (х, у)) = 0 (с) 9 (р, х). 


3. Подконтур, образуемый в С, цепочкой дуг (, Иа 


па’ 
о с оценками (м 9, ны 0, (уьу, }= 0, г= И р, и дугой 
(5), Си (рис. 2, г) устраняется заменой в матрице С’ действительного зна- 


чения Сс, Нах. Замена дает матрицу С? и список 0 (с ( у)) = 0 (с 9 (у, и). 


4. Чтобы устранить в 2, подконтуры, порождаемые цепочкой дуг [а 


а Изм] с оценками (и), у(у,. Уи) ®, 


а 0 2<г<5-1, 5+1 < г<р,и дугами (ри), ел (у, ), С ОЕ, 


Ск, С, #0 (рис. 2, д), образуем матрицу С, положив м, На бо 


матрице С’, 0 (с (х, у) = О (с 9 (ри) (уу). 
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Предположим, п, = | «|-2. Тогда с является значением оценки каждой дуги 
контура 7, за исключением (х, У) и (у, т 157. -Шх им (рис. 2, е). 
Если искомое решение содержит [9 У), то оно не содержит дугу (у, у, к иначе 
п, = [Е а|-1. Отсюда следует, что оценкой Ф ((х, У) с) является стоимость переста- 


НовКи Об (х, У) ‚ в которую входит элемент я У) и нет элемента (у, ый: к Другими 


словами, о, (5 У) представляет собой решение ЗН для матрицы С ”. полученной 
ГГ 
из С’ присвоением с элементу (уни) и элементам (х, у (1, У), уе. 6 


действительными значениями. 


Перестановка [©] [1% У) строится как и перестановка од [1% У), если по- 


Угу +1 
ложить а=у,, Р=У,.|, (а, Ь) — элемент, соответствующий (у, у, м Таким образом, 
при и, =[2.|-2 

2 
С (с и [о У)) + н, если построена перестановка © „, (5 У); 


фи У 19) = 
( } ) со в противном случае. 


Вершина ((х, у)°) не имеет допустимого продолжения, когда ф((х, у)°) = 
или построена циклическая перестановка с, (х, у). Если построена нециклическая 
перестановка с, (х, у), то вершине ((х, У) °) ставится в соответствие о, (х, У), 
 ((х, у) = © [а (*.›)). число контуров К (с,,(х,у)) в цикловом разложении 
перестановки с„, (х, У) и списки В (о, (1, у)}= В (6) (ху), 0 (въ (х,у))= 0 (в )ч 
® (а, ь). 

Таким образом, множество всех решений $ задачи рассматриваемого класса 
образует объедение непересекаю щихся подмножеств (+, ›)з) и ((х,у)). 


Всли обе перестановки бис (х, у) =. (х, у), построенные соответственно 


и * 
в вершинах ((. у) И ((х. у) °), оказываются циклическими, то $ =о0,, при ф (с - 


х 
ЗФ (с (3 у)) ив =50 [е у) при ф [6 (х, у)) < ф (с ыы Если получено только одно 
решение ЗН о,, или с [4% у), и оно представлено циклической перестановкой, то 
такая перестановка является искомым обходом. 
Концевые вершины, имею щие допустимые продолжения в строящемся дереве 
перебора, образуют множество активных вершин, которым соответствуют нецикли- 


ческие перестановки с оценками, равными стоимостям этих перестановок. Ветвлению ПОД- 
лежит активная вершина с наименьшей оценкой. Оценка в вершинеь, ДЛЯ которой Не 


* 
существует решения ЗН, равна ®. Построение обхода ® завершается нахождением 


концевой вершины и соответствующей ей циклической перестановки стоимостью, 
не большей, чем оценка в любой другой концевой вершине. 
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Рисунок 2 — Дуга с оценкой с изображена утолщенной линией, а дуга, исклю- 
ченная из контура 7, ,- пунктирной; а) для нециклической перестановки не 


существует соответствующего ей обхода, который содержит (х, У); исключается 
подконтур 6) (х, у, а) в) (^ НУ УНИ,.- У, ›--6 ух, г) (^ =Ур, У=У, 1 М, У, 


нЕ д) устраняются подконтуры (х У, УРУ,, М, М. Уз х) ‚ е) исключается 


контур Их. 


Модифицированный метод Литтла 


Перейдем к пошаговому представлению модификации метода Литтла. 


50. Из мат ИЦЫ = с- определить дважды приведенные матрицы Сис. 
р |, Пр р р ги С 


и соответствующие им суммы || и ц> констант приведения; 
если п! < ц> , то положить С =С,, ц=н , иначе С=С,, и=ц.; 
выполнить алгоритм решения ЗН для матрицы С; 
если для матрицы С, содержащей элемент с; =®, 1 ], ЗН не имеет решения, 


то конец: не существует обходов о 


если решением ЗН является циклическая перестановка о ‚ то - обход 8 
минимальной СТОИМОСТИ, 


с - нециклическая перестановка стоимостью фо =С(о)+ в в вершине ветвления 
Х (с ) ‚ представленной корнем дерева перебора, К (с) — число контуров в цикловом 
разложении перестановки о ‚ О(с<)=@, С.„„=С; т=1, 6„=о , Ф(о„)=С(о)+ы, 


К(с„)=К(с), 0(с„)=0(с), х={Х (с„}}. 


(соттет ‚ выполнение основного этапа начинается с ветвления корня) 
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51. Х=Х =[х (<„}} — список концевых вершин; процедурой ГАМС выполнить 


поиск в матрице жи элемента (^„, Ут) ‚ ИНИЦИирующего ветвление вершины Х (с т) , 


если выходом процедуры является 7 (2 Ут) = ©, то нет решений <, как со- 
держащих, так и не содержащих элемент (9х Ут ); перейти к шагу 55; 

сохранить результаты работы процедуры ГАМС: а Ут) 2 00, список 
в (с„)= ((, Гр (к,Г)= о), матрицу С’, множество всех элементов (Ё,Г), соответст- 


вующих дугам контура 7,, в котором содержится дуга (= м ), а также оценку 7 (К, Г) 


и 


каждой дуги (К, Ге А 21 ЕД, Д - цикловое разложение перестановки с | 


(соттеп!: формирование множества обходов, не содержащих [ее Ут )) 


+ 
Е Ко 


А = {х (с, }}; дая матрицы С', полученной из С’ заменой от 


Ут 


т Ут 


на © ‚ выполнить процедурой 5/\ построение перестановки сб, , ‚не включающей 
тит 


у 
элемент (хи, у }; 


если не построена перестановка о, ,‚ то вершина Х (с а ). порожденная 


Ут 


вершиной Х (с„), не имеет допустимого продолжения и получает оценку 


ф (с у ]= © ; перейти к шагу 53. Если построена циклическая перестановка с, , , 
т.т тит 


то Х [с ; — концевая вершина с допустимым решением =0 ‚ | СТОИМОСТЬЮ 
т Ут Хт Ут 


ф я )= С' (оз, )+ и ; перейти к шагу $3; 


— нециклическая перестановка, Х (с ) — активная вершина, которой 


Лт Ут Ат Ут 


Пе 
соответствует о, нижняя граница ф (с де )= С (с ии )+ и ‚ число контуров К с т 
в цикловом разложении перестановки о,,, (с, )=К(с„), 0(°„„)=0(с„) 5 


бы), 


53. Восстановить матрицу С’ в результате замены с, еС', на 
т 


= 00 
ти а 


действительное значение с 


АтУт Е 

(соттет : формирование множества обходов, включающих элемент (хи, Ут)) 

! 

54. Х=ХУ |х (с (ра )}}. Определить число и, элементов (К, 1) 2 с 

оценкой 7 (&,Г)=®; 
! у 

если и, =п,-1, то не существует перестановки с (хи, Уи ), включающей элемент 

[5 Уи (рис. 3, а); вершина х (с ее - )), порожденная вершиной Х (с > |? является 


Концевой с оценкой ф (с (х,,у, )) = ® ; Удалить матрицу С’ и перейти к 55. 


, о 
если п. < п. -2 ‚ ТО активной вершине Х (с [ее Ут )) Поставить в СООТВЕетствие 


перестановку с (ли, и}, нижнюю границу ф (о (хи, Уи )) = С'(6„)+н, С'(б„)= Сир (би), 
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К(с (хи. У„)) = К(о„), (с (х„,, у„)}=К(о„ Ух», и 0(с (хи, У) = О (9) Ру 


где содержимое Г, устанавливает один из четырех случаев исключения подконтуров 

ви. Ё=(унх,) фие: 3:6); 2) = (у, я | (рис. 3, в); 3) 6=(у„, м) (рис. 3, г); 

4) Г. = м и) (›,, т )ч (у, У) (рис. 3, д); удалить матрицу С’ и перейти к шагу $5. 

п, =, -2 (рис. 3, е); для матрицы ви полученной из С’ заменой на ® дейст- 
ели, (а,В)е 2 у (а, 6) = ®, ах 


т? 


вительных значений С.ь, С; ;, ]# Уи, С ни 


Ут’ 
р=у„, вызвать процедуру 5$ЛУ для построения перестановки о, (хи, Уи ), которая 
включает (хи, у) и не включает (4,6); 

если не существует 9, (хи, Уи), то Ф (с [А )) = < в концевой вершине 
Х (с [и Ут )); удалить матрицу С°; перейти к шагу $5; 


если сб, С У ) - циклическая перестановка, то Х (с (хи, Ут )) — концевая 


вершина, а 8 =су, (и Уи ) — обход стоимостью ф (с о Ут )) = С? (с г Е у ))+ 
+ ; удалить матрицу С". перейти к шагу 55; 

бар (и Уи ) — нециклическая перестановка; активной вершине Х [© [9 У )) 
поставить в соответствие перестановку бу, те Ут ), нижнюю границу ф (с [8 У ))= 


С’ (с р ))+ и, число контуров К (с а )) в цикловом разложении пере- 


СТановки ба в СПИСКИ В (с [5 Ут ))= в (си) О Уи }› 0( (х. у„))=0 („5 
(а, Ь); удалить матрицу С?. 


55. ш=ш+1; в списке Х найти вершину Х (<„) с наименьшим значением 


(с); 
* 
если © (о м 0 , ТО не существует искомого обхода © : конец; 


* 
если о;„ - циклическая перестановка, то 8 =Ох: конец, 


иначе сформировать матрицу С положив в исходной матрице С с, =®, 


сору 
11, су=® ,1#й, для всех (К, Г) Е (с) ис, = для всех элементов (1,5) 0 (с); 
перейти к шагу 51; 


“ * м 
Пример 1. Найти обход $ МИНиИмМальнои стоимости для матрицы 
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9 

со со со 
2 

со со со 
со со со со 


Выполним действия шага 50. Приведем матрицу С сначала по строкам, а затем 
полученный результат по столбцам, определив матрицу 


С — 


и сумму констант приведения щ =47+12+14+14+20+12+22+18+59+0+0+ 


+0+0+0+14+4=216. Приведем матрицу сначала по столбцам, а затем в получен- 
ной таблице выполним приведение по строкам. В результате получим матрицу 


© “миль 


и сумму констант приведения и, =61+12+14+14+20+12+51+18+10+0+0-0-6+0+ 
+0+3=221. Поскольку ц; < и›, положим С=С,, и=ц. Для матрицы С найдем ре- 
шение ЗН о = (1 7), (2,6), (3,8), (4, 5), (5,4), (6,2), (7,1), (8, 3), его стоимость С (о ) = ст + 
Соб + Сзв + С45 + 654 + 62 +сл1 + свз = 0+0+0+6+0+0-+0+6=6, и число К(С) контуров в 
циклическом разложении перестановки ©. Стоимость искомого обхода в’ ограни- 
чена снизу величиной фу =С (с) + ни = 6+216=222. Вершиной ветвления Х (С) яв- 
ляется корень дерева перебора. Сформируем матрицу С„„„, = С ‚ дублирующую на шаге 
50 матрицу С; т=1, 9, =6, ф(с,)=222, К (с, )=4, 0(°,)=0, Х = {Х (с, . 

51. Х=Х —{х (с ‚}} =@. Выполним процедуру ГАМС с матрицей Сор На 
входе: у (1,7) = 0+23-0=23, у (2,6) =23+27-0=50, у (3,8) =21+75-0=96, у (4,5) =23+ 
+12-6=29, у [5 4) =22+0-6=16, у (6, 2) = 13+21-0=34, у (7; 1) = 0+7-0=7, у (8, 3) =7+ 
+23-0=30. На выходе процедуры получим Г(3, 8) =96, К(с :) =©, С=С, 
2: ={(3,8), (8,3), п =2. 
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52. Х = {Хх (8). Положив с; = в матрице С’, получим матрицу С/, для 
которой процедурой 5 найдем решений ЗН оз = ((1,7), (2,6), (3,4), (4, 8), 
(5,2), (6,1), (7, 5), (8,3)), его стоимость С'(св)= =0+0+0+75+20+13+0+0=108 и 


нижнюю границу Ф(озв) = 108+216=324 в активной вершине Х (озз). Цикловое 


разложение перестановки ©.; состоит из двух контуров ((1, 7), (7,5), (5,2), (6), 
(6, 1), ((3, 4), (4, 8), (8, 3), К (озз ) = 2; К (оз) =9, О(сзв)= ((3,8)). 

$3. Восстановим матрицу С” заменой с; =< на сз; =0 в матрице С". 

$4. Х=Хы {х (с (3,8)}. Число элементов с оценкой © в 7; п; =0, щ=и;-2. 


| . . 
Положив в матрице С’ с; =®ис;, =®, ] #8, су =ю, 13, получим матрицу 


Процедура ЭЛ строит для матрицы С перестановку сз (3,8) = ((1,7), (2,6), 
(3,8), (4,3), (5,4), (6,2), (7,5), (8,1)), С” (свз (3,8))=7, х (свз (3,8)) =7+216=223. 
Цикловое разложение перестановки сб; (3,8) представлено контурами (1,7 | (17, 5), (5, 4), 
(.3). 6.8)... (2,6). 6,2), кб (,8))=2; в (6.8) = (8), 06 (5,8) = (83) 
удалим матрицу В 

95, =... В списке: Х = (х (6) Х (с (3,8))) вершина Х (с (3,8)) активна, 
имеет наименьшую границу ф (с (3, 3)) = С” (сз (3, 8) и = 223 и поэтому является 
вершиной ветвления; о) 558 Из исходной матрицы С образуем матрицу 


С ору, В КОТОрой С:=®, ] #8, СВЕ®, #3, Св=х. 
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$1. Х =(х (озз)) — список концевых вершин. После выполнения процедуры 
ГАМС получим 

у (1,7)=27, у(2,6) =50, у (3,8) =®, у (4,3) =46, у(5,4) = 63, у(6,2)=33, 
у (7,5) =6, у (8,1) = 

Г(5,4) =63, К(о›) = ((3,8), (8,1); 2, =((1,7), (7, 5), (5,4), (4,3), (3,8), (8, 

1)) ‚ ПП =б. Матрицы С’ на выходе процедуры ГАМС. отличается от С ору Тем, что в 
С’с=ю, 1. 

52. Х = (Х (сзз), Хх (954). В матрице С’ заменим с: =0 на с. Для получен- 
ной матрицы С' выполним процедурой 5ГУ построение перестановки ©.4, не 
содержащей (5,4). Процедура строит циклическую перестановку о = = ((1 т 
(52). (2,4, (4,6,(63,(38),(81) стоимостью ф(о54)=352. Вершина Х (054), поро- 
ждаемая вершиной Х (с (3, 8)) ‚ не имеет допустимого продолжения. 

$3. Заменив в матрице С" С54 =50 на 0, вновь получим матрицу С’. 

54. Х = (х (св), Х (и), Х (с (5,4))). п =2<п -2=4. Вершине Х (с (5.4)) 
соответствует перестановка с; (3,8), © (с (5,4) = 223, К (с (5,4))= К (сзз (3, 8))=2, 
В (с (5.4) = ®(вз (3,8) (5,4)= (3,8), (5,4), 0 (в (5.4) =0(в(38)) 545) = ((&, 3.45] 
(случай [| устранения подконтура в контуре 0, ). 

$5. ш =3. Ветвление выполняется в активной вершине Х (с (5; 4)) списка Х. 
Вершине Х (с (5.4)) соответствует наименьшая нижняя граница ф (с (5,4))=223; 


бз=0 (5, 4). Чтобы получить матрицу С для перестановки с}, в матрице С по- 


сору 


ЛОЖИМ СЗ;=®, ] #8, св =®, 13, су; =9, сц =®,{#5. Су, сея. 


(ов) а 
у (2,6)=50, 1(3,8)=®, а 46, 1(5,4)=ю, 1(6,2)=®, у 
6)=5 
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К(сз) = ((3, 8), (5, 4), (6, 2), (8, 1): 2. = (53 6), (6, 2) в цикловом разложении 
перестановки о (5,4) =озз (3,8), и› =2, п> =1. 

52. Х = (Х (сз), Х (5), Х (0). В С’ заменить съ =0 на о, чтобы получить 
матрицу С;. Результатом выполнения процедуры 5 для матрицы С! является ци- 
клическая перестановка о. (р т (7, 5), (5, 4), (4, 6, (6, 2). (2; ВНЕ 8), (8, 1)) стоимо- 
стью ф (05 ) = 273. 

53. Определим матрицу С’, заменив в © Со6 =°0 на 0. 

$4. Х = (х [бах (бя (с (2,6))}. Перестановка о (2,6) не сущест- 
вует, поскольку п, =и) -[. Вершина Х (с (2% 6) получает оценку ©. 

55. м =4. Вершине Х (с,,) соответствует циклическая перестановка с мини- 
мальным значением ф (с, )= 273. Поэтому искомый обход 8’ =.) СТОИМОСТЬЮ 
о (°)= С (еп }+ = 2735. 


Рис. 3. отображает процесс построения обхода 2’. П 


Ф(свз (3,8)) = 223 


Ф (свз (3,8) = 223 


Рисунок 3 - Дерево перебора для нахождения обхода в’=о.к; б| = ((1,7 

(2,6), (3,8), (4,5), (5,4), (6,2), (1,1), (8,3); озв = (1,7), (2,6), (8,4), 

(4,3), (5,2), (6,1), (7,5), (8,3); оз (3,8) = (1,7), (2,6), (3,8), (4,3), 

(5,4), (6,2), (1,5), (8,1); = (17), (7,5), (5,2), (2,4), (4,6), (6,3), (3,8), (8,1); 
бав = (17), (7, 5.6, 4), (4,6), (6,2), (2,3), (3,8), (8,1); (2,6) = 6 


) 
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5-3 А.Ю. Левченко, А.В. Морозов, А.В. Панишев 


Выводы 


Предложен модифицированный метод Литтла для решения задач класса комми- 
вояжера. Процедура ускоренного поиска нижней границы стоимости решения пред- 
ставляет собой вариант алгоритма ЗН, который строит ее решение из перестановки, 
полученной на предыдущем этапе ветвления. 
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А.Уи. Геусйепко, А.У. Мотгосотъ, А.У. Ратзйеу 
Месйатут ор Сотршаноп5 Ассеегайоп т [ле’; Метоа 
Гот эоГлиа Тгауейпе 5ще;тап С1а55 Та5К5 

Зошвоп$ зеагсН Гог пауейпе заезтап ргоет с1аз$ ш фе Ыпагу Бгапст® зсВете оЁ 
ЬгапсВ-апд-Боип тефо4 сап Ъе пойсеае ассеегаеа Бу теРеглие © а а аеогфт Юг 


зоуше а уапапё оЁР Фе азу1етитепе ргоМет (АР), изе {о сотрще 1о\ег Боип@$ Гог Фе 
Натотап топе” со$Ё [1]. Азяептепе ргоМеп?$ орйта! ретишайоп ог Фе тезШипе 


таблх сап Бе юип4 ш @те о ("?) ито зобоп Нот Фе ргеулои$ БгапсЬте $асе. 'ТВе Бгапсь 


ап Боип тефо4 фа и5ез а Газ а1гот т Гог зоуше Ше а$з1еттепе рго ет аз а 1юо\ег 
Боипа оРе зоаНоп?$ с05ё 15 ргорозе4. 


Статья поступила в редакцию 10.05.2012. 
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